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ABSTRACT 

A  two  dimensional  model  of  an  anchor  cable  subjected  to 
a  uniform  current,  fixed  at  the  sea  bed  and  moving  with  a 
specified  motion  at  the  surface  is  developed. 

A  steady  state  solution  is  developed  and  used  as  an 
initial  condition  for  the  numerical  solution  of  the  dynamic 
equations.   The  method  of  characteristics  is  used  to  carry 
out  the  numerical  integration. 
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TABLE  OF  SYMBOLS 

A  =  cross  sectional  area  (in2) 

An'  =  cable  acceleration  in  the  normal  direction  (ft/sec2) 

As '  =  cable  acceleration  in  the  tangential  direction  (ft/sec2) 

Cfi  =   drag  coefficient  for  a  cylinder 

d  =  diameter  of  cable  (in) 

D1  =  drag  force  per  unit  length  of  cable  (lbs/ft) 

Ys  '  =  depth  measured  from  the  mean  water  line  (ft) 

E  =  modulus  of  elasticity  (lbs/in2) 

K  =  cable  spring  constant  (lbs) 

L  =  length  of  cable  (ft) 

M  =  mass  per  unit  length  of  cable  (lb-sec2/f t2) 

M1  =  virtual  mass  per  unit  length  of  cable  (lb-sec2/f t2) 

p  =  density  of  seawater  (lb-sec2/f  f* ) 

V  =  angle  between  the  horizontal  and  a  tangent  to  the 

cable  at  a  point  on  the  cable 

S'  =  length  along  the  cable  measured  from  the  fixed  point 
on  the  sea  bed  (ft) 

T1  =  tension  (lbs) 

tQ  =  time  scale  (period  for  one  cycle)  (sec) 

TQ  =  steady  state  tension  of  the  cable  at  the  surface  (lbs) 

t1  =  time  (sec) 

V0  -  velocity  scale  (velocity  amplitude  at  the  surface) 
(ft/sec) 

V  =  current  velocity  (ft/sec) 

Vs '  =  tangential  cable  velocity  (ft/sec) 

V  '  =  normal  cable  velocity  (ft/sec) 


W   =  wieght  per  unit  length  of  cable  (lbs/ft) 

Ys'  =  vertical  distance  between  the  fixed  point  and  the 
top  of  the  cable  (ft) 
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I.   INTRODUCTION 

The  study  of  anchor  cable  dynamics  in  recent  years  has 
been  of  some  importance  to  those  interested  in  the  motion 
of  vessels  in  seaways.   The  mooring  of  buoys  for  navigation 
or  the  mooring  of  large  tanks  to  hold  oil  or  other  sea- 
transported  goods  are  further  examples  of  the  importance  of 
this  study.   The  motion  of  ships  riding  at  anchor  in  limited 
areas  would  be  of  interest,  for  example,  to  shipping  compa- 
nies that  service  small  ports. 

In  a  survey  [Ref.  l]  done  by  Casarella  and  Parsons,  it 
was  mentioned  that  the  first  experimental  work  on  mooring 
cables  was  probably  done  in  1917  by  Relf  and  Powell.   Since 
then  several  people  have  worked  on  the  mathematical  model 
of  this  problem. 

Through  the  years  the  model  has  changed  to  some  extent. 
Some  of  the  extensions  considered  have  been  the  inclusion 
of  tangential  drag  on  the  cable,  the  effect  of  surface  waves, 
the  effect  of  current,  the  effect  of  added  mass,  and  also 
the  effects  of  cross  currents  in  the  three  dimensional 
models  [Ref.  l].   Most  of  the  studies  done  have  been  carried 
out  for  the  static  configuration;  however,  dynamic  studies 
.with  arbitrary  motion  prescribed  at  the  surface  end  of  the 
cable  have  also  been  carried  out. 

With  the  advent  of  fast  computers,  time  required  to  get 
numerical  solutions  has  been  greatly  reduced;  however, 
computer  time  is  still  considerable  as  will  be  shown  in  this 
analysis . 

A. 
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Wilson  [Ref.  2]  arrived  at  a  model  very  similar  to  the 
one  used  in  this  thesis.   Wilson  considered  the  effect  of 
surface  waves  as  well  as  tangential  drag  on  the  cable; 
however,  the  computer  program  used  to  solve  the  problem  was 
proprietary  to  the  company  which  supported  his  research  and 
was  not  published. 

For  the  present  analysis,  the  drag  normal  to  the  anchor 
cable  only  is  considered  and  the  effects  of  waves  are 
neglected.   The  assumption  of  zero  tangential  drag  has  been 
included  in  the  general  approach  to  the  problem  and  seems 
to  be  valid  for  reasonable  currents  and  cable  configurations. 
Surface  waves  are  also  neglected  in  this  analysis  as  far  as 
drag  is  concerned.   The  current  applied  is  steady,  uniform 
and  in  the  horizontal  plane. 

In  the  dynamic  study,  the  method  of  characteristics  is 
used  to  solve  the  equations  of  motion  using  a  digital  computer, 


II.   THEORY 

The  problem  under  consideration  is  depicted  in  Figure  1. 
A  cable  of  diameter  d  and  weight  W  pounds  per  foot  is  anchored 
on  the  sea  floor  and  is  suspended  from  a  point  on  the  surface 
as  indicated  in  Figure  1.   A  uniform  current  of  magnitude  V 
is  directed  in  the  positive  X-direction  and  the  motion  of 
the  upper  end  of  the  cable  may  be  specified  in  some  suitable 
manner. 

Writing  Newton's  Law  for  an  elemental  length  of  cable 
as  depicted  in  the  insert  of  Figure  1,  the  following  dynamic 
equation  is  obtained: 

T'AY  -  WAS1  cos  V    -    D'As'   =   M'An'As'  (1) 

where  An'  is  the  normal  acceleration. 

Dividing  equation  (1)  by  As'  and  taking  the  limit  as  As'-*0 
yields : 

T1  1Z_   -  W  cos  T-D'=M'An'  (2) 

8s' 

The  term  M1  represents  the  virtual  mass  of  the  cable  per  unit 
length  and  D'  is  defined  as  the  cable  drag  per  unit  of  length. 
T*  denotes  the  tension  in  the  cable  and  the  angle  4'  denotes 
the  angle  between  the  tangent  at  a  point  on  the  cable  and  the 
horizontal.   In  summing  the  forces  in  the  tangential  direc- 
tion the  added  mass  is  taken  as  zero  and  as  a  result  the 
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FIGURE  1   DEFINITION  SKETCH 
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actual  mass  is  simply  the  cable  mass.   Also,  the  tangential 
drag  is  omitted  as  stated  in  the  introduction.   From  the 
forces  shown  in  the  insert  of  Figure  1,  the  equation  of 
motion  for  the  tangential  direction  takes  the  form: 

-T'  +  (T"  +  ill  As')cos  AY  -  WAs'sin(Y  +  £1)   =  MA  'As'  (3) 
9s '  2  s 

Again,  dividing  through  by  As1  and  taking  the  limit  as  As'+O 
yields: 

i^L  -  w  sin  V   =  MA  '  (4) 

3s'  s 

It  is  necessary  also  to  develop  expressions  for  the 
normal  and  tangential  accelerations  in  terms  of  the  normal 
and  tangential  velocities.   Referring  to  Figure  2,  the  normal 
and  tangential  acceleration  components  as  well  as  the  corre- 
sponding velocity  components  can  be  written  in  terms  of  the 
x  and  y  components  as 


As'   =   X  cos  V  +  Y  sin  V  (5) 

An'   =   Y  cos  Y  -  X  sin  *  (6) 

Vs'   =   Y  sin  *  +  X  cos  4"  (7) 

V  *   =   Y  cos  ¥  -  X  sin  *  (8) 

•      • 

in  which  X,  Y,  X,  Y,  denote  the  velocities  and  accelerations 
in  the  X  and  y  directions,  respectively.   Taking  the  partial 
derivatives  of  equations  (7)  and  (8)  with  respect  to  time 
and  using  the  definitions  of  V  '  and  V  '  as  given  in  equations 
(7)  and  (8)  it  is  possible  to  show  that  the  normal  and 
tangential  accelerations  may  be  written  as: 
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The  primes  indicate  that  these  terms  are  dimensional  terms. 
Dimensionless  forms  of  these  parameters  will  be  derived  later, 

Substituting  equations  (9)  and  (10)  for  the  acceleration 
into  equations  (2)  and  (4) ,  the  following  expressions  for  the 
two  equations  of  motion  are  obtained: 


T' 


3Y 
3T7 


-  W  cos  4"  -  D1 


M'  (IXlL.  +  vs 
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,  94<  ) 


(ID 


li-  -  W  sin  Y   =   M  (d  s   -  V 
3s1  3t' 


8Y 
3t' 


(12) 
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FIGURE  2   Coordinate  Systems 
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From  the  dynamic  equations,  (11)  and  (12) ,  it  can  be  seen 
that  there  are  four  dependent  variables,  V  '  ,  V„ ' ,  T'.V  ,  all 
of  which  are  functions  of  s'  and  t'.   In  order  to  solve  the 
equations  of  motion  two  other  equations  are  necessary.   These 
are  the  two  kinematic  relations  whose  derivations  follow. 

With  the  assumption  of  a  linear  relationship  between  the 
stress  and  strain  and  referring  to  Fig.  3,  it  is  possible  to 
write  the  following  relationship  between  the  rate  of  increase 
of  tension  and  the  rate  of  strain, 


Vn'  + 


4 

Ti 

+    3T« 
3s« 

As' 

3V     ' 
n 

,    As' 

~T 

7       > 

V  '  +  §«.  As' 

3s' 


FIGURE  3   Tangential  Forces  on  an  Element  of  Cable 
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which,  upon  taking  the  limit  as  As'->-0,  becomes 
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(14) 
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Equation  (14)  states  that  the  rate  of  change  of  tension 
equals  the  "spring  constant"  times  the  rate  of  stretching. 
In  the  case  of  a  solid  cable  of  uniform  cross  section,  K  is 
simply  AE  where  A  is  the  cross  sectional  area  and  E  denotes 
the  modulus  of  elasticity.   In  the  case  of,  for  example,  a 
chain  or  stranded  cable  the  value  of  K  would  become  somewhat 
more  complex  and  would  most  likely  be  determined  experimentally, 

v  ■  +  — —  as'    yv  •  +  L.  as' 

n    as'        /  s    3s«. 


FIGURE  4   Cable  Kinematics 

A  second  kinematic  relationship  may  also  be  introduced 

in  the  following  manner.   The  angular  velocity  of  a  segment 

As'  of  the  cable  may  be  specified  in  terms  of  the  velocity 

of  its  end  points.   Referring  to  Figure  4,  the  angular 

velocity  may  be  written  as: 

,  +  9VnJ_  As,  _    , 

il_  =  - 111 !L+Vs'iL_^i       (15) 

3t'  As'  9s1  As' 

which  becomes  upon  taking  the  limit  as  As'+O, 

3Y       3Vn'        9T 

=   +  Vs'  (16) 

3t'      3s'  3s" 


15 


The  four  equations  which  have  been  developed,  equations 
(11) ,  (12) ,  (14) ,  (16) ,  represent  the  complete  system  and 
are  all  that  are  required  to  determine  the  four  dependent 
variables  Vn',  V  ',  T' ,V    ,  in  terms  of  the  functions  s'  and 
t'.   In  summary,  the  basic  system  of  equations  is: 

T'  11—   -  W  cos  V  -  D1   =   M'  (2Xn_!_  +  V  '  11_)        (17) 
9s'  3t'      s   8t' 


111  -  W  sin  *   =   M(9Vs'  -  V  '  11—)  (18) 

8s.  at'     n  at' 

I  2ZL    =  £V  _  v  9j_  (19) 


K  9t'      9s'  8s 

9Y    _   8Vn'  .  ,,  ,  9^ 


+  Vs'  12—  (20) 

gt'      9s'         9s1 

The  four  equations  which  have  now  been  specified  equa- 
tions (17)  through  (20) ,  may  be  integrated  numerically  to 
obtain  a  solution.   It  is  possible  to  carry  out  this  numerical 
integration  using  a  finite  difference  formulation  in  the 
s'-t1  plane  with  a  rectangular  grid.   However,  it  is  also 
possible,  and  more  convenient,  to  use  the  method  of  charac- 
teristics as  outlined  by,  for  example,  Crandall  [Ref.  3], 
Using  his  suggested  method  it  can  be  shown  that  the  system 
is  hyperbolic  and  equations  (17)  through  (20)  can  be  trans- 
formed into  the  form  of  first  order  differential  equations 
valid  along  their  respective  characteristics. 

According  to  Crandall,  four  other  relationships  in 
addition  to  equations  (17)  through  (20)  are  needed  to  carry 
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out  the  integration  by  use  of  the  method  of  characteristics. 
These  relations  are  directional  derivatives  for  the  dependent 
variables  V_ ' ,  Vc * ,  T1  and  f  as  defined  by 


DVn'   = 


DVS' 


3  V  *        3V„  ' 

6    n   ds'  +    n   dt' 


3s' 
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Equations  (21)  through  (24)  denote  incremental  changes 
associated  with  small  displacements  ds '  and  dt'  in  the 
t'-s'  plane.   In  order  to  determine  the  characteristics  of 
the  set  of  equations,  equations  (17)  through  (24)  are  set 
in  a  matrix  form  as  follows: 


0 

M' 

0 

0 
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0 

3Vn'/3s' 
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-W  sin  Y 

-1 

0 

0 

0 
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0 
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0 

0 

0 

-K 

0 

Vn'K  0    0 

1 

3VS '/3t' 

0 

ds' 

dt' 

0 

0 

0    0    0 

0 

3V/3S1 

DV 

0 

0 

ds' 

dt' 

0    0    0 

0 

3*i'/3t' 

dvs' 

0 

0 

0 

0 

ds'  dt'   0 

0 

.3T'/3s' 

D41 

0 

0 

0 

0 

0    0   ds1 

dt' 

3T'/3t' 

DT' 

(25) 
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Expanding  the  determinant  of  the  square  matrix  and 
equating  it  to  zero  the  following  result  is  obtained: 


ds'4  -  (  -  +  — )ds,2dt'2  +  11*   dt'4   =   0  (26) 

M    M*  MM' 


Dividing  by  dt'4  equation  (26)  becomes: 


(dJLl)4-  (  *   +  11   )(^)2+  11*      =      0  (27) 

dt'       M    M'    df     MM' 


Equation  (27)  is  a  quartic,  the  solution  to  which  will 
produce  four  real  characteristic  directions. 
Letting 

Cs'2   =  -  (28) 

M 


Cn'2   =  11  (29) 

M" 


and  using  these  definitions  equation  (27)  becomes: 


(SSI)  -  (cs'2  +  cn'2)  (||4-)2  +  cs,2cn'2  =  0      (30) 


Solving  now  for  ds'/clt.'  gives  the  four  solutions: 
ds' 


=   ±Cn 
dt'       n 


(31) 


ds* 

—  =   ±CQ'  (32) 

dt"       s 
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These  slopes  define  the  four  characteristics  in  the  s'-t' 
plane  associated  with  the  system.   The  values  Cs'  and  Cn' 
agree  with  those  of  Wilson  [Ref.  2].   The  fact  that  there 
are  four  real  characteristics  classifies  the  system  as 
hyperbolic. 

Having  now  determined  the  slopes  of  the  characteristics, 
the  next  step  is  to  determine  the  first  order  ordinary 
differential  equations  along  each  of  these  characteristics. 
To  do  this,  the  right  hand  side  of  the  system  matrix 
equation  (Eq.  25)  is  substituted  into  any  column  of  the 
square  matrix  and  this  determinant  is  expanded  and  set 
equal  to  zero.   The  first  determinant  is  given  in  eq.  (33) . 
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Expanding  this  determinant  the  following  result  is  obtained: 
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where  C  ' ,  and  C  '  are  defined  by  equations  (28)  and  (29) 
and  C  =  ds' . 

at1 

It  may  be  noted  that  when  C  =  ±Cn'  both  sides  of  the  equation 
vanish.   When  C'  is  substituted  for  C  the  following  two 
equations,  valid  along  the  ±CS '  characteristic,  are  obtained: 


D(V')        D(¥)     1    D(T')      W  sin  T 

(+Cs")f +  V' -+ =   (35) 

Dt'  Dt'     MC  '   Dt1  M 


D(Vg')        D(f)     1    D(T')     -W  sin  V 

(-Cs'), Vn'  -+ =   (36) 

Dt1  Dt'     MCS'   Dt'  M 


Equations  (35)  and  (36)  are  the  ordinary  differential 
equations  valid  along  the  ±CS'  characteristics;  one  for  the 
positive  slope  and  one  for  the  negative  slope.   It  is 
interesting  to  note  here  that  the  Cs'  slopes  are  constant 
and  therefore  do  not  vary  throughout  the  complete  s'-t' 
plane. 

Equations  for  the  C  '  slopes  are  determined  in  a  similar 
manner  by  substituting  the  right  hand  side  of  the  matrix 
equation  system  (Eq.  (25))  into  another  position  in  the 
square  matrix;  the  fifth  column  will  produce  the  desired 
result.   After  expanding,  the  following  equation  is  obtained: 


(Cs'2-c2)   ^Li+c(C2-C  2)   E<2L+V  .(c   .2^2,   DJhO, 
Dt'  Dt'  S  Dt' 


=       (Cs'2-C2)(-   W   cos    f    4-    D'j        (3?) 

M' 
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In  this  case,  if  C  is  replaced  by  ±CS'  the  equation 
becomes  trivial.   However,  if  ±Cn'  is  substituted  for  C 
the  following  equations,  valid  along  the  second  character- 
istics Cn',  are  obtained: 


(+Cn)    ^^  (Vs'-  Cn')  ^r     «   -  (  W  COSM,       )   (38) 

(-Cn)    D(Vn'}  +  (VB'  +  CM  21.  =   -  (  W  cos  *  +  D'  )   (39) 
Dt'       s     n    Dt'  Mi 


Equations  (35)  ,  (36)  ,  (38) ,  and  (39)  represent  four 
ordinary  first  order  differential  equations  which  are  valid 
along  their  respective  characteristics.   The  total  deriva- 
tives - —  are  defined  according  to   D  .  =  jl ±C  8 


Dt'  Dt'    3t'     9S' 

ds' 
where   C  =  -rrr  denotes  the  slope  of  the  characteristic  in 

the  s'-t'  plane.   Since  these  equations  are  first  order 

equations  they  are  easily  integrated  along  the  characteristic 

directions. 

With  a  uniform  current  directed  as  shown  in  Figure  1, 

the  drag  per  unit  length  acting  normal  to  the  cable  may  be 

written  in  terms  of  the  fluid  velocity  relative  to  the  cable 

in  the  form 

D*   =   JspCdd(V  sin  V  +  Vn')|V  sin  t  +  Vn'|     (40) 

where  p  denotes  the  fluid  density,  C-,  denotes  the  drag 
coefficient,  V  denotes  the  current  magnitude  and  d  denotes 
the  cable  diameter  or  frontal  area  per  unit  length.   The 
absolute  value  is  used  in  order  to  maintain  the  sign  of  the 
relative  velocity. 
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A.   DIMENSIONLESS  PARAMETERS 

The  four  characteristic  equations  (35) ,  (36)  ,  (38) ,  and 
(39)  can  now  be  put  in  a  dimensionless  form.  Various  terms 
need  to  be  defined  before  the  actual  analysis  is  done. 

The  solution  to  the  problem  as  depicted  in  Figure  1 
depends  on  a  number  of  parameters.   That  is,  the  periodic 
motion  at  the  surface  may  be  characterized  by  a  period  (t  ) 
and  a  velocity  amplitude  (V0) .   Also,  at  steady  state,  the 
tension  at  the  surface  T  as  well  as  an  angle  between  the 
horizontal  and  a  tangent  f   to  the  cable  may  be  specified. 
The  unstretched  length  (L)  of  the  cable  can  be  used  to 
normalize  all  length  scales  and  the  velocity  amplitude  may 
be  used  to  normalize  the  velocities.   Applying  these 
definitions  to  the  characteristic  equations  the  following 
dimensionless  parameters  may  be  defined: 

(41) 


(42) 


t 

= 

f 
fco 

vs 

= 

V 

Vo 

Vn 

= 

v  ■ 
n 

vo 

T 

= 

T' 
To 

V 

= 

V 
Vo 

cn 

= 

<V 

L 

(43) 


(44) 


(45) 


(46) 
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C_   =   C  '  (JL)  (47) 

S       S    L 


(48) 


B    =   -2_2  (49) 


*t0 

Vo 

Vo 

LMVQ 

tow 

V0M' 

VoPCddto 

(50) 


D    = (51) 

2M' 


E    =   —  (52) 


(53) 


(54) 


Vo 

F 

= 

L 
vo*o 

R 

_ 

K 

WL 


,RA^ 
s 


Cc   =   (f^)  2  (55) 

F 


C    =  (*£L)h  (56) 

n      VAF 


Applying  the  dimensionless  groupings  (41)  through  (56) 
to  equations  (35)  ,  (36)  ,  (38)  ,  and  (39)  ,  the  characteristic 
equations  can  be  written  in  the  following  dimensionless  form: 


(+CS)    -  D(Vg)  +  Vn  E11L  +  B_  DJT)_   =   A  sin  ,           (57) 

Dt  Dt  Cs  Dt 

D<vs)  D(¥)  B  D(T) 

(-C_) Vn  i—-*"  +  —  ^T~ L      =    "A  sin  V           58 

S       Dt  n  Dt  cs  Dt 
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D(Vn)  DY  _ 

(+C)   —  +  (Vs-FCn)  —  =  -  C  cos  V  -  DLE  sin  ¥  +  Vn]  |E  sin  *  +  vj 

Dt  Dt 

(59) 


D(VJ  DV  r  -,. 

(-C^)   —  +  (V-+FCJ  —  =  -  C  cos  f  -  D[E  sin  ¥  +  V.]  E  sin  f  +  V 

n    Dt       s  ^1  Dt  n  n1 

(60) 


Equations  (57)  through  (60)  may  now  be  solved  numerically 
on  the  digital  computer  to  obtain  the  dynamic  solution. 
However,  in  order  to  proceed  with  the  solution  it  is  first 
necessary  to  determine  an  initial  condition  as  well  as 
appropriate  boundary  conditions  on  the  surface  and  ocean 
bottom. 

B.   STEADY  STATE  SOLUTION 

The  steady  state  solution  may  be  obtained  by  setting  Vn, 

Vs  and  _  equal  to  zero  in  equations  (57)  through  (60) . 
at 

Accordingly,  the  following  steady  state  form  of  equations 
(57)  and  (58)  is  obtained: 

B  *L  =  A  sin  f  (61) 

ds 

which  is  equivalent  to  equation  (12)  if  the  parameter  defini- 
tions are  substituted.   Equations  (59)  and  (60)  become: 

FCn2  $1     =   C  cos  ¥  +  DE2sin2  f  (62) 

ds 

Equation  (62)  is  equivalent  to  the  equation  of  motion  for 
the  normal  direction,  equation  (11) .   Since  the  characteris- 
tic equations  agree  with  the  equations  of  motion  as  far  as 
steady  state  is  concerned,  equations  (59)  and  (60)  will  be 
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the  governing  equations  for  the  steady  state  solution. 
Also,  the  fact  that  both  sets  of  equations  agree  is  an 
indication  that  the  characteristic  equations  were  derived 
correctly. 

Substituting  equation  (56)  for  C   and  dividing  equation 
(61)  by  equation  (62)  provides  the  following  expression: 

dT      sin  ¥  d4» 


T      cos  f  +  DE   sin   Y 


(63) 


This  equation  can  be  integrated  analytically  between  the 
limits  of  T  to  TQ  and  4*  to  4'  .   Upon  integration  the  follow- 
ing result  is  obtained: 

(2d  cos  4-0-lVl+i^)  (2DEi  cos  T  _i+Vl+f2^) 
2  2  V         C  C2     C  C2 

T  U+~^2";  =  — I /   2  4 " 7   2  4      (64) 

CSf.  cos  *0-l<VW&JL)  (Si  cos  *  -lVl-^4-) 
C      °       r2     C  C2 


If  there  should  be  zero  current  then  equation  (63)  would 
reduce  considerably  since  the  parameter  E  would  drop  out. 
An  analysis  is  done  in  Appendix  B  to  compare  equation  (64) 
without  the  current  term  to  the  solution  obtained  directly 
from  equation  (63)  setting  E  equal  to  zero. 

Once  T  is  known  from  equation  (64)  then  the  following 
expression  can  be  used  to  get  s: 


ds   =   BCT  dy (65) 

AC  cos  4"  +  ADE2  sin2  1 


25 


Equation  (65)  can  not  be  analytically  evaluated,  and  there- 
fore it  was  integrated  numerically  on  the  digital  computer. 

C.   BOUNDARY  CONDITIONS 

With  the  steady  state  solution  used  for  the  initial 
values  only  boundary  conditions  at  the  surface  and  sea  bed 
need  be  specified.   At  the  sea  bed  the  normal  and  tangential 
velocities  will  be  zero  since  the  cable  is  anchored.   With 
these  two  relations  plus  equations  (58)  and  (60) ,  the  four 
dependent  variables  can  be  determined  along  the  sea  bed. 

At  the  surface,  equations  (57)  and  (59)  can  be  used; 
however,  two  other  equations  need  to  be  specified  such  that 
the  four  dependent  variables  can  be  calculated.   These  two 
relations  can  be  specified  in  several  ways.   For  example, 
the  surface  velocities  may  be  written  as: 

Vs   =  Vy  sin  ¥  +  Vx  cos  ¥  (66) 

Vn   =  Vy  cos  ¥  -  Vx  sin  T  (67) 

where  V"x  and  V  may  be  specified  functions  of  time. 

The  two  relationships  required  at. the  free  surface  may 
also  take  the  form  of  the  two-dimensional  equations  of 
motion  for  a  floating  vessel.   In  this  case  tension  would 
be  related  to  the  position,  velocity  and  acceleration  of 
the  floating  object. 
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III.       NUMERICAL    PROCEDURE 

Equations    (57)    through    (60)    represent   four   ordinary   first 
order   differential   equations  valid   along   their   respective 
characteristics.      In  order   to  carry  out   the   numerical   integra- 
tion,   these   equations   are   placed   in   the   following   difference 
form: 

-Vs5+Vnl¥5+(CL)T5    =  Asin*lDt-Vsl+VnlV<§r)Tl  (68) 

s         .  s 


VnC.+  (V^-FCr0)T[.     =      (-C  cos  y--D[E  sin  Y9+Vn9]|E  sin  Y-+VLJ  )Dt 


'n5    lvsi2  rv"n2'T5 


2   vn2- 


2  vn2' 


Wri2+(Vs2-pCn2),1'2 


.(69) 


Vh_+(V,,,+PC:T1,)Yc     =     (~C  cos  T--D[E  sin  t_+Vn,]|E  sin  H'_+Vm|)Dt 


'n5'  vvs3  iWn3' '5 


3     "3- 


3     n3' 


wn3+fv«'wcn3)¥3 


Yks-^'S+^^S  =  "A  Sin  *4  Dt+Vs4-vn4,F4+(i-)T4 
s  s 


(70) 
(71) 


Figure  5   System  of  Characteristics 
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The  subscripts  1,2,3,4  denote  the  initial  point  where  it 
is  supposed  that  all  quantities  V  ,  V  ,  T,  and  4*  are  known. 
In  order  to  "march  out"  and  determine  new  values  at  the  new 
location  denoted  by  the  subscript  5,  the  four  linear  equations 
(68)  through  (71)  are  solved  simultaneously.   Letting  the 
right  hand  sides  of  equations  (68)  through  (71)  be  equal  to 
AA,  CC,  DD,  and  BB,  respectively,  the  resulting  expressions 
for  the  dependent  variables  at  the  new  location  5  are: 

T5   =   (CC  -  DD/((VS2  -  Vs3)  -  F(Cn2  +  Cn3))        (72) 
T5   =   (AA  +  BB  +  4-§(Vn4  -  Vnl))/(g£)  (73) 

V     =   BB  +  V   w   -  (E_ )t  (74) 

VS5     aa   r  114*5    ^c  '  5  K      ' 

s 

Vn5   =   CC  -  (VS2  -  FCn2)^5  (75) 

The  basic  method  of  calculation  indicated  in  equations 
(72)  through  (75)  was  used  to  determine  values  of  the  depen- 
dent variables  V  ,  V  ,  T,  and  *f  throughout  the  t-s  plane. 
The  numerical  scheme  utilized  consisted  of  first  dividing  the 
s  axis  into  equal  subdivisions.   The  solution  for  V  ,  V  ,  T, 
and  V    at  t=0,  or  initial  condition,  was  obtained  from  the 
steady  state  solution  of  equations  (64)  and  (65) .   The  steady 
state  solution  provided  values  for  V  ,  V  ,  T,  and  T  at  all  of 
the  nodal  points  along  the  S  axis.   The  layout  in  the  t-s 
plane  is  shown  in  Figure  6. 

Using  equations  (72)  through  (75)  a  new  set  of  values 
were  then  determined  at  the  next  row  of  points  located  at 
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time  Dt  =  Ds/C   from  the  s  axis.   Once  all  the  values  of  the 
s 

four  dependent  variables  were  determined  at  this  value  of 
time,  the  solution  was  marched  out  another  time  step,  Dt,  so 
that  values  of  the  dependent  variables  were  eventually 
obtained  for  a  large  area  on  the  t-s  plane. 

It  may  be  noted  in  Figures  5  and  6  that  the  ±C   charac- 
teristics do  not  originate  from  a  nodal  point  where  values 
of  the  dependent  variables  are  calculated.   In  the  numerical 
procedure  the  method  of  calculation  involved  a  linear  inter- 
polation between  points  1  and  4  in  Figure  5  in  order  to 
evaluate  the  dependent  variables  at  points  2  and  3. 

As  indicated  previously,  the  steady  state  solution  was 
utilized  to  determine  the  initial  values  of  the  four  dependent 
variables  along  the  s  axis.   On  the  t  axis,  corresponding  to 
the  ocean  floor,  it  is  also  necessary  to  apply  some  kind  of 
boundary  condition.   In  the  example  worked  out,  this  boundary 

condition  was  taken  as  V  =  V  =0  which  corresponds  to  a 

n    s  r 

fixed  anchor  point. 
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FIGURE  6   Layout  of  Characteristics 
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IV.   DISCUSSION  OF  RESULTS 

In  order  to  work  a  numerical  example,  it  is  necessary  to 
specify  values  for  the  parameters  appearing  in  equations  (41) 
through  (56) .   For  the  example  worked,  a  steel  cable  of  uniform 
cross  section  with  a  diameter  (d)  of  2  inches  and  a  modulus 

of  elasticity  (E)  of  30  x  10   psi  was  selected.   The  density 

3 
of  the  steel  was  480  lbs/ft   and  the  cable  was  1000  feet  long. 

The  period  (t  )  of  the  motion  at  the  surface  was  10  seconds 

and  the  velocity  amplitude  V  was  3.0  feet/second.   The 

2    4 
density  of  seawater  was  taken  to  be  2.0  lbs-sec  /ft  .   The 

Trd2 
virtual  mass  as  described  by  M'=M  +  C  p  —. —  where  C   denotes 

J  m    4  m 

the  added  mass  coefficient  [Ref.  2],  and  p  the  density  of 
seawater,  was  calculated  to  be  3.49.   The  added  mass  coeffi- 
cient used  in  the  above  calculation  was  taken  as  1.5.   The 

drag  coefficient  C,  was  taken  as  1.5  which,  according  to 

5 
Wilson  [Ref.  2],  is  valid  for  Reynolds  numbers  less  than  10  . 

It  is  now  possible  to  calculate  values  for  the  parameters  in 
equations  (48)  through  (56)  . 

With  an  accurate  steady  state  solution  for  tension  having 
been  calculated  by  equation  (64)  ,  a  check  on  the  numerical 
integrations  done  in  the  characteristics  method  is  possible. 
That  is,  it  is  possible  to  apply  steady  state  boundary  condi- 
tions at  the  free  surface  and  recompute  the  steady  state 
solution  using  the  method  of  characteristic  procedure.   These 
results  can  then  be  compared  with  the  original  steady  state 
solution. 
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In  order  to  do  this  the  steady  state  boundary  condition  at 
the  surface  was  written  in  the  following  manner: 

T  cos  *|    =cos  ¥  (76) 

s=l       o 

This  equation  states  that  the  horizontal  component  of  tension 
at  the  surface  was  held  constant.  The  vertical  component  was 
written  as: 


T  sin  ¥    =sin  <F   +  C  (Y   -  Y   )  (77) 

1  s=l       o    os    s  ' 

Equation  (77)  is  a  simplified  expression  representing  the 

vertical  component  of  tension  as  caused  by  a  massless  vessel 

floating  on  a  still  water  surface  in  water  of  depth,  Y 

o 

The  dimensionless  water  depth  Y    is  defined  as  Y    = 
,  *  s  s 

\  o  o 

/"  sin  ¥ds  I., n   and  likewise  the  vertical  span  of  the  cable 

0  •  U  1 

at  time  t  is  defined  as   y   =  /   sin  ¥  ds | .  .   The  product 

s  ■  0  t 

C  (yc  -  y   )  represents  a  buoyant  force  caused  by  a  sinking 

o 

of  the  vessel  a  distance  of  (y   -  y   ) .   The  constant  C   is 

o 

used  to  describe  the  interaction  of  the  ship  with  the  free 

surface.   The  "spring  rate"  C   represents  the  slope  of  the 

buoyant  force-sinkage  curve  and  is  accordingly  defined  as 

C   =  (A  , ) Y   ) L/T    where  A  ,  is  the  water  line  area,  v 
o     wl   sw  '    o         wl  sw 

is  the  specific  weight  of  seawater,  and  L  is  the  length  of 
the  cable.   Using  equations  (76)  and  (77)  as  the  surface 
boundary  condition,  the  method  of  characteristics  was  "marched 
out"  several  time  steps  in  order  to  compare  the  results 
obtained  from  the  method  of  characteristics  procedure  with  the 
initial  steady  state  results. 
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The  results  of  the  above  described  procedure  are  shown 
in  Figure  7.   The  results  corresponding  to  t=0  were  calculated 
by  use  of  equations  (64)  and  (65)  .   Using  these  initial  values, 
the  method  of  characteristics  procedure  was  then  used  to 
continue  this  steady  state  solution  for  some  time.   It  may  be 
seen  from  Figure  7  that  some  drift  in  the  results  took  place 
due  to  the  numerical  integration  by  the  method  of  characteris- 
tics when  compared  to  the  initial  solution.   The  figure  shows 
that  the  vessel  sank  slightly  as  the  tension  increased.   This 
change  from  the  initial  solution  is  not  to  be  construed  as 
inaccuracy,  however.   The  new  steady  state  solution  simply 
represents  a  second  possible  steady  state  solution. 

Since  the  steady  state  solution  using  the  method  of 
characteristics  apparently  provided  good  results,  it  was 
decided  to  add  a  dynamic  boundary  condition  at  the  surface  in 
order  to  induce  a  dynamic  response  of  the  cable.   The  boundary 
conditions  at  the  surface,  equations  (76)  and  (77) ,  were 
therefore  modified  slightly  to  include  motion.   Equation  (76) 
was  applied  so  that  the  horizontal  component  of  tension  was 
held  constant.   The  vertical  component  of  tension,  however, 
was  expressed  in  the  following  way: 

T  sin  V  I   .   =   sin  ¥   +  C  [  (Y  -Y   )  +  A   sin  2irt]   (78) 
1 s=l  o    oL   s   s      m 

o 

Figure  8  depicts  the  physical  problem  defined  by  the 

boundary  condition  equation  (78) .   The  equation  describes  the 

vertical  component  of  tension  as  caused  by  a  massless  vessel 

floating  in  a  wave  of  elevation  A   sin  2iit  measured  with 

3  m 
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FIGURE  8   Dynamic  System 
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respect  to  the  mean  water  level.   The  term  A  represents  the 
wave  amplitude  and  is  made  dimensionless  with  the  cable  length, 
The  terms  within  the  square  brackets  represent  the  displace- 
ment of  the  vessel  from  its  mean  position  with  respect  to  the 
free  water  surface.   As  defined  previously,  C   denotes  the 
slope  of  the  buoyant  force  versus  sinkage  curve. 

Figure  9  is  a  graph  of  the  dynamic  response  of  the  cable 
to  the  above  described  motion.   On  this  graph  are  plotted  the 
dimensionless  tension  and  depth  against  a  dimensionless  time 

scale.   Tension  is  made  dimensionless  with  T  which  is 

o 

defined  as  T   =  WL/sin  V      ,  ¥   =  45°  being  the  initial  steady 
state  angle  at  the  surface.   The  depth  is  made  dimensionless 
with  the  cable  length,  which  is  1000  feet  for  the  example 
worked.   Finally,  the  time  is  made  dimensionless  with  the 
period  of  the  motion  at  the  surface  which  is  10  seconds.   The 
T   calculated  from  the  above  expression  is  approximately 
14,150  pounds.   The  results  of  the  steady  state  analysis, 
Figure  7  is  also  shown  on  Figure  9.   Figure  9  shows  a  sinu- 
soidal response  which  is  expected  from  the  boundary  condition 
driving  the  motion.   Although  the  average  tension  at  the 
surface  should  be  near  1.0,  it  can  be  noticed  that  after 
about  three  cycles  the  tension  appears  to  be  oscillating 
around  a  mean  value  which  has  a  6  per  cent  shift.   This  new 
mean  position  is  equal  to  the  one  calculated  from  the  method 
of  characteristics  under  steady  state  conditions. 

Finally,  it  may  be  noted  from  Figure  9  that  a  phase  shift 
exists  between  the  position  of  the  upper  end  of  the  cable  and 
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the  tension.   Specifically,  the  tension  tends  to  lag  the 
displacement  somewhat.   This  lag  is  caused  by  the  dynamics 
of  the  phenomena;  a  juxtaposition  of  static  solutions  would 
necessarily  show  the  tension  and  displacement  to  be  in  phase. 
The  amount  of  computer  time  required  to  obtain  numerical 
results  depends  on  two  factors,  the  number  of  subdivisions 
along  the  cable  and  the  slope  of  the  C   characteristic.   In 
general,  the  number  of  calculations  required,  and  consequent- 
ly, the  computer  time,  is  directly  proportional  to  the  square 
of  the  number  of  subdivisions  along  the  s-  axis  or  the  cable 
length.   The  size  of  the  time  step  advanced  in  each  computa- 
tion is  related  to  the  slope  of  the  C   characteristic 

c  s 

according  to  Dt  =  Ds/C   where  Ds  is  fixed  by  selection  of 

the  number  of  subdivisions  of  the  cable  and  C   is  defined  as 

s 

C   =  (K/M)   .   Thus,  in  order  to  reduce  the  required  computer 
time  it  is  necessary  to  use  the  minimum  number  of  subdivisions 
of  the  cable  that  will  give  accurate  results  and  use  small 
values  of  C  .   Unfortunately  C   is  dependent  on  the  cable 
physical  properties  and,  accordingly,  is  fixed. 

In  view  of  the  definition  of  C   it  is  apparent  that  the 
time  steps  for  very  "stiff"  cables  with  large  elastic  con- 
stants K  will  be  small.   However,  for  lighter  cables  with 
smaller  elastic  constants  such  as  would  correspond  to  nylon 
rope,  the  time  steps  will  be  larger  and  the  calculation 
procedure  more  efficient.   For  the  particular  example  worked, 
30  subdivisions  on  the  cable  were  used  and  it  took  approxi- 
mately 30  minutes  of  computer  time  to  compute  2.9  cycles 
corresponding  to  29  seconds  in  real  time. 
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I, 
The  value  of  C   =  (K/M)  2  has  the  physical  significance 

of  being  the  speed  of  propagation  of  a  stress  wave  along  the 
cable  while  the  second  characteristic  C   =  (T'/M1)  2   has 
the  physical  significance  of  being  the  speed  of  propagation 
of  a  transverse  wave  along  the  cable.   Accordingly,  for  all 
physically  realistic  cables  the  stress  wave  speed  would 
normally  be  much  larger  than  the  transverse  wave  speed  and 
consequently,  C  would  be  greater  than  C  .   In  the  case  of 
steel  cables,  the  stress  wave  speed  is  normally  quite  large 
making  the  time  step  small.   Moreover,  in  the  limiting  case 
of  the  inextensible  cable  the  calculation  procedure  breaks 
down  because  the  time  step  becomes  zero.   However,  in  this 
case  or  in  cases  where  the  cable  stretching  can  be  disregarded, 
the  problem  actually  becomes  much  simpler.   That  is,  the 
system  of  four  equations  (57-60)  are  replaced  by  the  two 
equations  (59)  and  (60)  which  are  valid  along  the  ±C   charac- 
teristics only.   The  ±C   characteristics  are  lost  completely 
and  the  system  reduces  to  a  much  simpler  two-characteristic 
system. 


39 


V.   CONCLUSIONS 

The  mathematical  model  presented  in  this  thesis  provides 
a  method  for  the  solution  of  cable  dynamic  problems.   Although 
a  simplified  boundary  condition  was  employed  on  the  surface 
for  purposes  of  generating  numerical  results,  the  results 
obtained  indicate  that  the  method  is  valid  and  can  be  used  in 
connection  with  any  type  of  boundary  condition  on  the  surface. 

Excessive  computer  time  is  the  greatest  problem  to  be 

overcome  in  producing  numerical  results  and  two  methods  for 

reducing  these  were  discussed.   In  cases  where  stretching  can 

be  disregarded  the  computer  time  can  be  reduced  by  orders  of 

magnitude  by  utilizing  the  two-characteristic  model  since 

the  time  step  is  in  this  case  defined  by  Dt  =  Ds/C   and  C 
c  J  n       n 

is  normally  small  compared  to  C  .   In  cases  where  the  cable 
is  relatively  stiff  and  yet  long  enough  such  that  stretching 
cannot  be  disregarded,  no  solution  to  the  problem  of  excessive 
computer  time  presents  itself. 
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RECOMMENDATIONS 

The  following  recommendations  are  advanced  on  the  basis 
of  experience  gained  during  the  course  of  the  present  study: 

1.  A  surface  boundary  condition  should  be  applied  that 
contains  all  the  terms  involved  in  a  surface  vessel's 
motion  including  mass,  added  mass,  damping  and  buoyant 
force. 

2.  The  two-characteristic  model  should  be  developed  and 
compared  to  the  solutions  using  the  four-characteristic 
method  in  order  to  evaluate  the  effect  of  stretching. 

3.  A  study  of  the  effect  of  cable  subdivision  size  on 
accuracy  should  be  carried  out. 

4.  The  model  appearing  in  this  thesis  should  be  extended  to 
include  tangential  drag,  and  surface  wave  drag  on  the 
cable  and  should  be  extended  to  the  three  dimensional 
case. 
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APPENDIX  B 

It  is  interesting  to  note  that  if  the  current  is  zero 
then  equation  (63)  becomes: 

dT   _   sin  V  d$  ,__, 

T~  "   cos  v  (78) 

This  equation  can  be  easily  integrated  into  the  following  form: 

cos  V 
T   = -  (79) 

COS  Y 

By  setting  the  current  equal  to  zero  in  equation  (64)  it  is 
not  apparent  that  the  result  will  be  the  same  as  equation  (79) . 
It  is  possible  to  show  that  equation  (64)  does  reduce  to 
equation  (79)  by  the  following  method. 

4  2 

T  4-4--  4E*DZ 

Letting      e  =   5— 

cr 

and  using  a  binomial  expansion  for   (1+e)    the  following 
is  obtained: 

(1+eP   =  1  +  ^e   -  I  e2  +  Jg.  e3  (80) 

Substituting   equation    (80)    into  the   numerator   of   equation    (64) 
produces   the   following   numerator: 


to        1J5  m  12  1       3W      1      h  „    .    1  12.13. 

(-    2    -  ^    cos   YQ   -    ^    e      -   ^   £    ) (-    2    ^    cos    T    +  f    "    g    E      +   T6    £    ) 

(81) 
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The  coefficient  of  the  e  2  term  will  produce  cos  V    since  the 
other  terms  will  be  of  higher  order.   A  similar  analysis  can 
be  done  in  the  denominator  producing  a  coefficient  of  cos  ¥  . 
Equation  (64)  now  reduces  to: 


cos  ¥ 

T   =   ^  (82) 

cos  ¥  v   ' 


which  is  identical  to  equation  (79) . 
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